Highly nonlinear dynamics of third-harmonic generation by focused beams 
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Beams that experience third-harmonic generation (THG) also experience Kerr effects. With 
C**) i Kerr effects, beams do not take simple Gaussian shapes, but exhibit nonlinear dynamics. These 

. nonlinear dynamics have an effect on the THG accumulated by focusing and then diverging beams. 

We formulate a self-consistent and complete set of nonlinear Schrodinger equations for a pair of 
coupled beams - a fundamental and its third-harmonic. Numerical simulations show that the Kerr 
^ . nonlinearities allow some third-harmonic to propagate to the far-field even for zero or negative 

Q ■ phase-mismatch. This is because the nonlinear dynamics break the beams' reflection symmetry 

•"^ ' about the focal plane and therefore increases far-field THG by changing some of the interference 

from destructive to constructive. THG conversion efficiencies are computed as functions of several 
■ beam parameters. 
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PACS numbers: 42.25.Ja, 42.65.Jx, 42.65.Sf, 42.65.Ky 
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I. INTRODUCTION 
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^ . We study the dynamics of two-color beams in a nonlinear isotropic medium. We take the initial conditions to 
^ ' be a one-color beam with Gaussian profile. With a cubic (x^) nonlinearity, which is the lowest order possible 
. , in an isotropic medium, the possible nonlinearities are THG, Kerr, and Raman. Any other nonlinearity requires 
^ ■ either a different susceptibility (second-harmonics, high-harmonics), or more than two slowly- varying envelopes (sum 
^h ^| and difference frequency generation). This is one of the simplest nonlinear opticsproblems, and is for that reason 
f-^ i important; it has an application to third-harmonic generation (THG) microscopy 

The usual model of THG takes the fundamental beam to be Gaussian, and (inconsistently, or as an approximation) 
has THG as the sole nonlinearity 0, El ■ This model has exact analytic solution in which the third-harmonic (TH) 
beam takes a Gaussian profile. For a phase-mismatch that is zero or negative, the energy in the TH peaks at the 
, focus of the fundamental beam, and, after the focus, destructive interference causes all the TH to be re-absorbed by 
■ the fundamental. The latter is often said to be due to the Guoy shift Q, the phase-shift of tt that a Gaussian beam 
experiences in going from a far-field, through a focus, to a far-field. But, in fact, the details show it to depend on 
the beam's phase and shape everywhere along its path. However, the fast (electronic) nonlinear response of optical 
materials that yields THG also generates Kerr effects (self- and cross-phase modulation). Therefore, one should 
include Kerr and perhaps also Raman effects when modeling THG. These can have quite drastic effects on a beam £|. 
Although it was recognized as early as 1973 that Kerr effects could influence THG j5|, and although Kerr, THG, and 
& dispersion have been examined in studies without transverse spatial dynamics 0, Q , the results of transverse nonlinear 
. 3h ' dynamics on THG have not heretofore been studied in quantitative detail. 

We have derived a set of coupled nonlinear Schrodinger (NLS) equations for two slowly-varying envelopes, with the 
relevant nonlinearities treated rigorously and consistently |g. Because of the microscopy application, we are interested 
in very tightly focused beams and pulses, with large momentum and frequency spreads. The NLS equations arc thus 
• • . given to all orders in dispersion and diffraction; numerical simulations are carried out with a method that is accurate 
. 1 to all orders relevant for the grid 0,0,01- We express the field as a fundamental, A w (x,t), centered about a carrier 
wave at frequency ujq, and a TH, As u (x, t), centered about a carrier wave at frequency Sloq, 
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In the NLS equations QJ, = d 2 + d 2 is the transverse Laplacian, n(ui) is the index of refraction at frequency w, 
and f3(ui) = n(w)w/c is the wavenumber. The first summation terms on the right-hand sides (RHS) of Eqs. Q are 
dispersion, the second summation is diffraction, the third summation is cross-dispersion/ diffraction. The nonlinear 
polarization includes THG, Kerr, and Raman effects 0,0,0. Self-steepening terms 0,0], of the form d{\A\ 2 A) / dt 

and self- frequency shifting terms of the form ^S—A are contained in the first-order time derivatives of the nonlinear 
polarizations P u (x,t) and P 3uj (x,t) in Eqs. (JIJ. The NLS equations JTJ show that self-steepening terms and self- 
frequency shifting terms are just the first terms of a family of higher order nonlinear terms. 
The nonlinear polarization is taken to be of the form 

P u [x,t) = 3x clcc (-^o;-^o,-^o,3^ ) exp{-i[3/?M-/3(3^o)]2}^(x,t)* 2 ^(x,t) 

+ 3x clcc (-^o;^o,-wo,wo) |Au(x,i)| 2 A.(x,i) +6x olcc (-w ;3wo,-3wo,wo) \A 3u (x, <)| 2 A,(x, t) 

r OC 

+ A,(x,i) / [3x nucl (-^o;^o,-^o,^ ; s) \A u (x,t- s) | 2 + 3x nucl (-^o;^o, -3w ,3wo; s) \ A 3u) (x,t-s) \ 2 }ds 
Jo 

/>oo 

+ A 3uJ (x,t) exp(-2iw s)3x nucl (-w ; 3^o, -^o, -3w ; a)A U) (x,t - s)A 3uJ (x,t - s)*ds (2a) 
Jo 

P 3w (x,t) - x clcc (-3w ;w ,wo,^o) exp{i[3/3M-/3(3w )]z}A w (x,i) 3 

+ 6x cloc (-3w ;a;o,-wo,3wo) |A w (x, i)| 2 A 3w (x, i) + 3x olcc (-3w ; 3uj q , -3uj q , 3lu ) \A 3uJ (x,t)\ 2 A 3u (x,t) 



A 3uj (x,t) / [3x nucl (-3w ;w ,-w ,3a;o; s) 14, (x, t - s)\ 2 + 3x nucl (-3w ; 3w , -3w , 3w ; a) |A 3w (x, i - s)\ 2 ]ds 
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P OC 

A^(x,t) / exp(2w s)3x nucl (-3w ;^o,-wo,3wo; s) ^(x,i-s)M 3w (x,i-s)ds . (2b) 



This breaks up the nonlinear response into an electronic (fast) part and a nuclear (slow) part. It also assumes that 
the electronic part of the response may be considered frequency-independent on the scale of the pulse bandwidths; 
it does not make this assumption for the nuclear part of the response. For the calculations, we take a more specific 
nonlinear polarization, essentially the standard model of the third-order susceptibility of fused silica [l^, Ha, Ua, 
plus a generalization, 

P w (x, t) = exp{-i[3/3(w ) - 0(3<J O )}z} x THG 4,(x, i)* 2 A 3w (x, t) + x° lcc (|^(x, t)\ 2 + 2\A 3lu (x, i)| 2 )A,(x, t) 
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+A u (x,t) / X nucl (s)(\AUx,t-s)\ 2 + \A 3LJ (x,t~s)\ 2 )ds 
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+A 3LJ (x,t) / cxp(-2iw s)x nucl (s)^(x,t - s)A 3LU {x,t - s)*ds , (3a) 
Jo 

P 3a) (x,t) = exp{i[3^(wo) - /3(3w )]«} ix™ G ^(x, i) 3 + x° loc (2|^(x, t)| 2 + |A 3w (x, t)| 2 )A 3w (x, t) 
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+A 3w (x,t) / x nucl (s)(|^(x,t- S )| 2 + |A 3w (x,t- S )| 2 )d S 
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+A u (x,t) / cxp(2 lWo .s)x nucl (s)^(x,t-.s)*A 3w (x,t-.s)d.s. (3b) 
Jo 

Here the THG coefficient is decoupled from the other electronic susceptibilities; this is outside the usual model 
for pure fused silica in which the electronic contribution is considered instantaneous and the nuclear contrubution 
takes the form arising from a single damped harmonic oscillator, x(t'i ti, t 2 , £3) = x clcc ^(^ — ti)6(ti — t 2 )6(t 2 — t 3 ) + 
X nucl (^i -t 2 )5{t-tx)5{t 2 -t 3 ), with x olcc = n 2 (l-/ Raman ), x nucl (0 = n 2 f Raman (T 2 + t 2 )t^t 2 2 exp(-i/T 2 ) sin(t/n), 
n 2 = n(w)cn^/(2?r), n\ = 2.8 x 1CT 20 m 2 /W, /Raman = 0.18, n = 12.2 fs, and r 2 = 32.0 fs. If the electronic THG 
susceptibility differs from the other electronic susceptibilities, the electronic contribution to the susceptibility is not 
instantaneous compared to all scales. Direct experimental measurements of the THG susceptibilities are available 
[l8L IT9II20I I . The more recent measurement of the THG coefficient ^ THG is smaller than the electronic contribution 
to self-phase modulation x° lcc by a factor of almost four; older measurements [2(j give a THG coefficient smaller 
than the electronic part of the self-phase modulation by a factor of about 1.5. In the absense of direct experimental 
measurements of all the nonlinear polarization coefficients for doped silica, we use the simplest case x T = x° lec 
in numerical simulations, but discuss how the results scale for different values of the THG coefficient. Vector effects 
are neglected. Inter-band Raman scattering is negligible at the carrier frequencies because of the fast relative phase 
oscillation; but in the pulse simulations, inter-band Raman scattering is possible between the lower frequencies within 
the higher frequency band, and the higher frequencies within the lower frequency band. 



3 



II. NUMERICAL SIMULATIONS 



A numerical NLS propagation scheme may be said to be accurate to all orders of dispersion and diffraction if it 
is accurate to as many orders as there are grid points. Accuracy up to all available orders requires the index of 
refraction over the entire numerically represented frequency range. In the split-step fast Fourier transform scheme, 
linear propagation is carried out in momentum space; the algorithm may be made accurate to all orders by putting 
the frequency dependent index of refraction directly into the formulas for the propagators ||, i.e., 



>(3w + w) 2 - k\ - 0{3u> o ) 



HM*,m k^- { V '^'} (4a) 



Here T and T 1 are Fourier and inverse Fourier transforms in cc,y,i, k± = y£; 2 + fc 2 is the transverse momentum, 

and uj is frequency. The linear dispersion effects are contained in (3{ui) = n(uj)uj/c; diffraction is contained in k±. The 
offsets by loq and 3wo are due to the fact that the fields are slowly-varying envelopes about carrier waves at those 
frequencies; when these appear, uj is the frequency relative to the offset. There is no need to compute the coefficients 
for dispersion, diffraction, self-steepening, etc., explicitly because they are contained implicitly in the linear dispersion 
P{u>). Since we simulate focusing and collapse, during which diffractive, dispersive, and nonlinear length scales can 
easily change by factors of a hundred or more |l2l [l4j , we allow the propagation step to vary such that it remains an 
order of magnitude less that the smallest relevent scale. 

The effects that we are interested in are larger and clearer when the fundamental and TH are not too far from 
being phase-matched. Ak = 3/3(wo) — (3(3ujq) ~ 0, or n(ujo) w n(3ajo). Optical materials may be doped to obtain 
desired properties [12, W$ ■ We consider silica doped with neodynium to obtain approximate phase-matching between 
the fundamental A — 1.5 and TH A = 0.5 fim. To model the frequency dependent index of refraction, we use the 
Sellmeier relation for fused silica j^, but add one additional resonance at ANd = 0.59 /im, which is neodynium's largest 
resonance in the vicinity of our TH. For this, phase-matching is achieved at the Sellmeier coefficient £?Nd ~ 0.0138, 
which corresponds to a few percent doping of the material. Phase-mismatch is varied by changing the dopant 
concentration (in this model, the Sellmeier coefficient i?Nd)- We take the conventional nonlinear coefficients for pure 
silica 1 2j , as given above. 

We simulated the propagation of both pulses and continuous beams, over a range of light intensities, phase-mismatch 
values, and focusing strengths. The effects of the nonlincarity are clearly visible in Fig.^ which shows the evolution 
of beams with varying intensity. For all simulations, we took initial conditions with zero power in the TH, and the 
fundamental in the form of a Gaussian with some radial phase factor Ai(x, t; z = 0) = A% cxp[— (ic\ + 2 In 2/Wi)(x 2 + 
y 2 )}. The radial phase factor c\ can be related to the wavefront radius of curvature R\ via the magnitude of the 
wavevector (3(u) = n(uj)uo/c according to the formula R\ = 0.5/3(w)/ci. Because the fundamental focuses very 
intensely, while only a relatively small part of it is converted to TH, we illustrate the peak intensity of the fundamental 
beam and the power of the TH beam. In Fig. the input beam has initial width FWHM = 50 /im, a radial phase 
factor c\ = 5.0 x 10 9 /m 2 that brings the beam to a focus in 0.6 mm, and the material is phase-matched (the reciprocal 
of the phase-mismatch is more than an order of magnitude greater than the simulation distance). The fact that the 
normalized TH power varies with fundamental input intensity shows that the beams experience nonlinear dynamics, 
changing shape as the intensity changes. In Fig. ^ the initial wave front radii of curvature are small and their effect 
on the position of the focus overshadow the intensity dependence. One can see that the beams approach the linear 
limit at low intensities, as these curves start to almost overlap. Fig. [21 shows the normalized far-field TH power as a 
function of input power, for the data in Fig. ^ and for another set of runs with weaker focusing. 

For a qualitative explanation of these figures, we first contrast our results with the model in which 
THG is the only nonlinearity 0, Q- Here, in the case with phase- matching, the TH power Pth(z) = 

(3/2)(27r/nc) 2 (u;/c) 4 (7rW r 1 2 ) 3 x^ 3 ^ I^IW + ( z ~ Zfocus) 2 /^ 2 ] reaches a maximum at the focus of the fundamental and 
then drops off as the inverse square of z/zi, where z\ = ttW 2 /Ao is the Rayleigh range of the fundamental and the 
TH beam, and W\ is the width of the fundamental at its focus. Our simulations show that at high intensities, the 
Kerr effect causes the beams to lose reflection symmetry about the focal plane. With this asymmetry, THG from the 
incoming and outgoing beams does not fully interfere destructively, and allows some TH to propagate to the far- field. 
The lower-intensity curves in Fig. ^ do not quite drop off as the inverse square of distance because the numerical 
simulations did not start out from minus infinity, but began with merely a large (finite) beam width. The reflection 
asymmetry in Fig.^is due partly to the nonlincarities and partly to starting with a finite initial beam width. As long 
as the amount of energy in the TH is relatively small, the nonlinear dynamic effects will remain in the fundamental 
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beam, and the TH peak intensities and beam powers may scale up or down by a uniform factor, but will be otherwise 
unaffected. Clearly, there is significant variation with the fundamental input power. Moreover, this variation is quite 
different for different focusing conditions. 

In another series of simulations, phase-mismatch is varied. Fig. [21 shows the peak intensity of the fundamental and 
the TH beam power with initial FWHM = 50 fixa, radial phase-factor c\ = 5.0 x 10 9 /m 2 , and phase-mismatch from 
-12.9 to 51.8 mm -1 . Because the TH power first increases and then decreases with phase-mismatch, we show the 
former range on one plot and the latter on another. Clearly, there is a residual far-field TH and its power depends on 
phase- mismatch. Fig. 0] shows the far-field TH conversion efficiency as a function of phase mismatch, for the runs in 
Fig.[3]and another set of runs with weaker focusing, c\ = 2.5 x 10 9 /m 2 . As in the linear model Q, THG by a focusing 
beam is maximized around a certain phase-mismatch; but the nonlinear dynamics complicate the results considerably. 

In a further series of simulations, the intensity and phase-mismatch were held constant, and the radial phase 
factor varied. This mixes up a few physical effects, since initial conditions with small radial phase factors cannot be 
considered as starting from the far-field. Fig. [5] has extremely small phase-mismatch (Afc = 0.067 mm -1 ), and radial 
phase factors which vary from zero (i.e., starting at a focus) up to ci = 5.0 x 10 9 /m 2 . The horizontal axis is on a 
log scale to help visually distinguish the superimposed sumulations. The peak fundamental intensity varies, but over 
the range studied, the peak TH power is rather insensitive to the wavefront radius of curvature of the fundamental, 
but the far-field TH power does vary considerably. Fig. [B] shows the far-field TH power as a function of radial phase 
factor, for the runs in Fig. and also for a series of runs with significant positive and negative phase-mismatch. In 
increasing (decreasing) the tightness of the focusing, the nonlinear and diffraction lengths both decrease (increase); 
the complexity of the results for far-field THG reflect the complexity of the nonlinear dynamics of the beam. 

Fig. H shows the dependence of the third-harmonic power on the when the THG suseptibility ^ THG is varied. The 
curve labeled ^ THG = ^°i ec corresponds to the same conditions used in Fig.[3]but with Afc = 0.067 mm -1 (i.e., almost 
phase- matched) . The third harmonic power scales with |x TH | 2 when the third-harmonic intensity is small since then 
the nonlinear dynamics of the fundamental is unaffected by TH and the generation of TH field is proportional to 
X THG . The values of x THG used in the calculations shown in Fig. 0are a factor of 1.0, 1/3.9, and 1/1.5 times the 
value used in the previous figures, where the latter two factors correspond to the measured values of x THG reported 
in Refs. OH III 

Numerical simulations showed the dynamics of pulses to differ from those of continuous waves in essentially two 
ways. First, the group velocities of the fundamental and TH will, except for special cases, not be the same. A TH 
pulse will thus generally walk off from a fundamental pulse. This tends to reduce destructive interference; it also limits 
the effective distance over which the pulses interact. Compared to beams, TH pulses tend to carry off a larger part of 
the fundamental pulse energy. Secondly, the dynamics of pulses (3+1 -dimensional) arc qualitatively different to those 
of continuous beams (2+1-dimensional). For strongly focused but not very short pulses, these differences tend to be 
minor. A thorough analysis of focused pulses in this system is quite involved, and will be pursued elsewhere. Thus, 
the detailed results for continuous- wave beams apply to pulses for which group- velocity differences are relatively small 
or for pulses that are relatively long. 

III. CONCLUSIONS 

In a nonlinear medium, intense beams or pulses of finite diameter which converge to a focus and then diverge may 
exhibit nonlinear dynamics that significanly affect propagation dynamics. These nonlinear effects break the beam's 
reflection symmetry about the focal plane. The greater the intensity, the bigger the difference between incoming 
and outgoing beams. THG with such an input beam or pulse produces TH in the far-field when the phase-mismatch 
between the fundamental and its third- harmonic is zero or negative (and small). We have quantitatively demonstrated 
this for several cases. When phase-mismatch is positive, where some far-field TH power is possible in the essentially 
linear case, the nonlinear beam dynamics complicate the accumulation of TH power. The non-zero far-field THG 
for zero or negative phase-mismatch - a qualitatively new effect for homogeneous media - is affected by, and thus 
contains information about the medium in the region of the beam focus. For an inhomogeneous medium, and for 
THG microscopy, these effects should be understood, either to be utilized or better avoided. 
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List of Figure Captions 

Fig. 1. Peak intensity of the fundamental (Ai = 1.5 /mi) beam, and power of the third-harmonic (A3 = 0.5 fxm) 
normalized by the cube of the input power. The initial conditions have a range of intensities (1.01, 4.05, 16.2, 36.4, 
64.7, 101, 122 and 145 KW, corresponding to circled points in the appropriate curve in Fig. [2J, but identical beam 
width, FWHM = 50 /mi, radial phase factor, c\ = 5.0 x 10 9 /m 2 (i.e., transverse phase exp[— ici(a- 2 + y 2 )]), and all 
are phase-matched. 

Fig. 2. far-field third-harmonic beam power, normalized by the cube of the input power, vs. input power. The 
medium has TH phase-matching. Two curves are shown, one the far-field results from Fig. 1, and another with 
smaller radial phase factor, that brings the beam to a focus in about 1.2 mm. 

Fig. 3 Peak intensity of the fundamental beam, and power of the third-harmonic. The curves represent beams 
with identical power, width, and radial phase factor, while the phase-mismatch goes from -12.9 to 51.8 mm -1 
(corresponding to circled points in the appropriate curve in Fig.^J. The middle plot shows the TH increasing with 
phase-mismatch, and the bottom plot shows TH decreasing with phase-mismatch. 

Fig. 4. Far-field conversion efficiency vs. phase-mismatch. Two curves are shown, one the far-field results from Fig. 3, 
and another with smaller radial phase factor, that brings the beam to a focus in about 1.2 mm. Initial intensities are 
held constant. 

Fig. 5. Peak intensity of the fundamental beam, and power of the third-harmonic versus propagation distance for 
different fundamental radial curvatures. The curves show beams with phase-matching and initially identical power 
and width; the initial radial phase factors go from zero (starting at a focus) up to 5.0 x 10 9 /m 2 , by increments of 
0.5 x 10 9 /m 2 (the curves here correspond to circled points in the appropriate curve in Fig. |SJ). 

Fig. 6. Far-field conversion efficiency vs. radial phase factor. Two curves are shown, one the far-field results from 
Fig. 5 (phase-matching), and another in a medium with positive phase- mismatching. 

Fig. 7. Peak fundamental intensity and third-harmonic power versus position z. for various values of x T ■ The third 
harmonic power scales with |x T | 2 when the third-harmonic intensity is small. 
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FIG. 1: Peak intensity of the fundamental (Ai = 1.5 pm) beam, and power of the third-harmonic (A3 = 0.5 (im) normalized by 
the cube of the input power. The initial conditions have a range of intensities (1.01, 4.05, 16.2, 36.4, 64.7, 101, 122 and 145 
KW, corresponding to circled points in the appropriate curve in Fig. but identical beam width, FWHM = 50 /»m, radial 
phase factor, ci = 5.0 x 10 9 /m 2 (i.e., transverse phase exp[— ici(x 2 + y 2 )]), and all are phase-matched. 



8 




10 3 10 4 10 5 

input power (in fundamental beam) [W] 

FIG. 2: Far-field third-harmonic beam power, normalized by the cube of the input power, vs. input power. The medium has 
TH phase-matching. Two curves are shown, one the far-field results from Fig.0 and another with smaller radial phase factor, 
that brings the beam to a focus in about 1.2 mm. 
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FIG. 3: Peak intensity of the fundamental beam, and power of the third-harmonic. The curves represent beams with identical 
power, width, and radial phase factor, while the phase-mismatch goes from -12.9 to 51.8 mm -1 (the specific values correspond 
to circled points in the appropriate curve in Fig. 0J. The middle plot shows the TH increasing with phase-mismatch, and the 
bottom plot shows TH decreasing with phase-mismatch. 




FIG. 4: Far-field conversion efficiency vs. phase-mismatch. Two curves are shown, one the far-field results from Fig. [3] and 
another with smaller radial phase factor, that brings the beam to a focus in about 1.2 mm. Initial intensities are held constant. 
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FIG. 5: Peak intensity of the fundamental beam, and power of the third-harmonic versus propagation distance for different 
fundamental radial curvatures. The curves show beams with phase-matching and initially identical power and width; the initial 
radial phase factors go from zero (starting at a focus) up to 5.0 x 10 9 /m 2 , by increments of 0.5 x 10 9 /m 2 (the curves here 
correspond to circled points in the appropriate curve in Fig.|^J. 




FIG. 6: Far-field conversion efficiency vs. radial phase factor. Three curves are shown, one corresponding to the far-field results 
of Fig. (phase- matching) , and two others corresponding to positive and negative phase-mismatch. 
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FIG. 7: Peak fundamental intensity and third-harmonic power versus position z, for various values of \ . The third harmonic 
power scales with |x THG | 2 when the third-harmonic intensity is small. 



